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We consider a tiieoretical model for the chiral smectic A twisted ribbons observed in assemblies 
of fd viruses condensed by depletion forces. The depletion interaction is modeled by an edge energy 
assumed to be proportional to the depletant polymer in solution. Our model is based on the Helfrich 
energy for surface bending and the de Gennes model of chiral smectic A liquid crystals with twist 
penetration at the edge. We consider two variants of this model, one with the conventional Helfrich 
Gaussian curvature term, and a second with saddle-splay energy. A mean field analysis of both 
models yields a first -order phase transition between ribbons and semi-infinite flat membranes as 
the edge energy is varied. The phase transition line and tilt angle profile are found to be nearly 
identical for the two models; the pitch of the ribbon, however, does show some differences. Our 
model yields good qualitative agreement with experimental observations if the sign of the Gaussian 
curvature or saddle-splay modulus is chosen to favor negative Gaussian curvature. 
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I. INTRODUCTION 

Smectic A (Sm-A) layers expel twist and bend defor- 
mations, just as magnetic fields are expelled from bulk 
superconductors [ll, • In analogy with the London pen- 
etration length, the typical distance to which a magnetic 
field penetrates into a superconductor, two penetration 
depths can be defined in Sm-A liquid crystals describ- 
ing respectively the penetration of twist and bend defor- 
mations at the edges or around isolated defects of Sm- 
A layers. When a Sm-A sample is composed of chiral 
molecules (molecules without mirror symmetry), twist 
deformations appear intrinsically, driving the formation 
of novel structures, due to the competition of the twist 
deformations with the tendency of the molecules to build 
a perfect twist-free smectic layer Q- When the layer 
forming tendency dominates, the twist is restricted to a 
band at the layer edge, but when the twist deformations 
are strong enough, simple flat membranes are replaced 
by a variety of twisted structures, including twisted rib- 
bons, double helices, arrays of twist walls in membranes, 
and periodic arrays of pores. In this paper, we study a 
theoretical model for one result of this basic competition 
of ordering forces, the transition between flat disks and 
twisted ribbons. 

Recently, there have been both experimental and the- 
oretical studies to visualize and quantify the effects of 
twist and bend deformations on chiral smectic A (Sm- 
A*) single-layer membranes in the weak chirality limit, 
as well as investigating the relative stability of different 
geometries of membranes seen in experiments. These 
membranes are formed under certain conditions in aque- 
ous solutions of filamentous virus particles and a non- 
binding polymer, which acts to condense the virus par- 
ticles into dense phases by depiction forces. Barry et. 
al. have studied flat Sm-A* monolayer membranes in the 
form of disks composed of rodlike fd virus particles to 



measure the twist penetration length A* at the edges of 
the disks [j|. The persistence length of these viruses is 
2.8 ± 0.7 ^m and the length of fd viruses was around 1 
/iTO, making them nearly rigid rods (5|. Because they act 
like hard rods, the condensed phases they form due to the 
depletion force are mainly entropy driven [6|, |7| ■ These 
viruses form cholesteric and Sm-A* phases with increas- 
ing concentration. The radii of the disks studied were 
tens of micrometers. Due to the length scale defined by 
the particle length, the twist penetration could be visu- 
alized using optical microscopy. Due to the tilting of the 
virus particles relative to the layer normal, in the twist 
penetration region, the twist could be measured by the 
resulting change in retardance of transmitted polarized 
light. The twist penetration length At was determined to 
be about 0.5 nm. Because the disks were relatively large 
compared to this, the twist penetration profile fit well 
to the analytic theory of twist penetration for a semi- 
infinite membrane, with the director tilting tangential to 
the edge of the membrane. From the estimated maximum 
tilt at the disk edge, the product of penetration length 
and cholesteric wave vector q was estimated as approxi- 
mately q\t — 0.71, for the particular sample studied. 

In analyzing the energy of the edge of a disk, there are 
two important terms, the bare edge energy, or line ten- 
sion, of the edge, 7, and the net reduction of free energy 
in the material near the edge due to the twist penetra- 
tion. The line tension is controlled by the concentration 
of depletant polymer in solution. As this concentration 
is lowered, 7 is reduced. In a simple analysis, if the sec- 
ond term then becomes larger than the first, then the net 
energy associated with the edge is negative, and the ma- 
terial will adopt a structure that maximizes the length 
of edge, relative to the area of fiat disk. This can be 
achieved in several ways. 

One mechanism, previously studied, is a transition 
from semi -infinite membranes to small disks as the struc- 



ture that minimizes the mean free energy density of the 
system. The mean free energy density of small Sm-A* 
disks, which have both twist and bend deformations at 
the curved disk edge, has been calculated Ql- It is found 
that small disks with a twist/bend penetration length on 
the order of the disk radius are at least metastable rela- 
tive to large disks upon reducing the line tension of the 
edge. The magnitude of the critical line tension at which 
a second-order phase transition occurs from semi-infinite 
membranes to small disks is affected by the magnitudes of 
the twist wave vector q and the twist penetration length 
A*. 

Rather than simple small disks, a diversity of struc- 
tures is observed experimentally to replace large disks 
as the concentration of depletion agents is lowered Q- 
Twisted ribbons (minimal surfaces to a double helix) are 
commonly seen. In this paper we carry out a mean-field 
study of the twisted ribbon structure. Throughout this 
work we assume that the smectic order parameter coher- 
ence length ^ is less than At, so that the Sm-A* phase 
is analogous to a type II superconductor. Additionally, 
we focus here on a model in which twist and bend de- 
formations near the edges of the Sm-A* layers are as- 
sumed to be driven by the chirality of the molecules, as 
a continuation of the previous work [J, Q . We obtain the 
twist penetration profile as a function of the width of the 
twisted ribbon. Furthermore, we are able to study the 
stability of these structures and determine the first-order 
phase transition from isolated Sm-A* twisted ribbons to 
semi-infinite flat layers. 

The present study is organized as follows: In the next 
section we introduce our model for the elasticity of smec- 
tic A* membranes and then apply it to twisted ribbons. 
In Sec. Ill, we use our model to study the transition 
between infinite disks and twisted ribbons. We offer con- 
cluding remarks in the final section of the paper. 



II. FREE ENERGY OF Sm-A* MEMBRANES 

A. Elasticity of membranes 

We model Sm-A* membranes using the Helfrich 
model [10, \lM for the surface bending energy and the 
de Gennes model [l|, 01 for the Sm-A phase generalized 
to include chirality. Our model is an extension of one 
introduced in Ref. |l2l| to include variations in the tilt 
angel of the director relative to the normal to the sur- 
face of the membrane. We write the free energy F of an 
Sm-A* membrane as follows: 



F ■ 



ifH + fn)dA + jIdl, 



(1) 



where fn and /„ are the Helfrich and de Gennes free en- 
ergy densities respectively, dA is the area element of the 
membrane, and dl is the arc element length of the edge. 
The last term represents an edge energy with energy per 



unit length (line tension) 7. The Helfrich energy density 
fn is given by: 



fH = \k{2Hf + ~kKG, 



(2) 



where H and Kq are the mean and the Gaussian curva- 
tures of the surface respectively, k is the bending rigidity 
and k is the Gaussian curvature modulus. There is no 
spontaneous curvature in Sm-A* membranes given their 
up-down symmetry, and thus we have not included such 
a term. Furthermore, twisted ribbons have zero mean 
curvature so henceforth we set i? = 0. In terms of N 
the local unit normal vector to the surface the Gaussian 
curvature is given by |13j : 



KG^]rV ■ (n(V • N) - (N • V)N 



(3) 



If the molecular director field n is everywhere parallel to 
the layer normal N, then the Gaussian curvature is equal 
to the Frank saddle-splay energy density [13, [ij]. In 
Sec. Ill Cl we consider an alternate model of the membrane 
where we replace the Gaussian curvature by saddle splay; 
we compare the results of the two models of a twisted 
ribbon in Sec. IIIII 

The de Gennes free energy density for the nematic- 
Sm-A* phase transition is given by: 
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Here ^ is the smectic order parameter, Ki,K2 and K^ 
denote the splay, twist, and bend Frank elastic constants 
respectively, q is magnitude of the spontaneous twist 
wave vector arising from molecular chirality, and 5\\ is the 
deviation of the director from the layer normal. In a bulk 
Sm-A phase twist and bend distortions which are incom- 
patible with constant interlayer spacing will be expelled; 
however, as noted by de Gennes both twist and bend can 
be introduced at layer edges or around defects by allow- 
ing local tilting of n relative to the layer normal. Two 
penetration lengths can be introduced, A2 = {K2/DY/'^ 
and A3 = {K^/D)^/'^, describing respectively the pene- 
tration lengths of twist and bend into the smectic phase. 
In the limit where the smectic order parameter ^ is 
uniform and nonzero (i.e., deep in the smectic A* phase) 
the de Gennes free energy density reduces (up to a con- 
stant) to: 



/„ = iifi (V • n)' + ii^2(n -Vxn-qf 



-Ksin X V X n)2 



2 



(5) 



where 9 is the tilt angle of the director with respect to the 
layer normal and the tilt energy density C is proportional 
to the coupling D appearing in Eq. (|3|). 



In the single Frank elastic constant approximation, 
K = Ki — K2 — K^ , which we will henceforth employ in 
the analysis of Eq. ([T|), /„ is given by: 
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(V • n)^ - 2qn • (V x n) + (V x n)^ + q^ 



+ -Csin^i 




(6) 



Mathematically, a membrane can be represented as a 
surface with position vector Y = Y(ui, U2) embedded in 
a three-dimensional Euclidean space, where ui and U2 
are real variables parameterizing the surface. To calcu- 
late the free energy of the membrane we need the follow- 
ing geometric quantities [l^ [ig : 



FIG. 1; The local membrane basis showing the tangent vec- 
tors Yi and Y2 and the normal vector N perpendicular to 
them. The director n is tilted by an angle 9 with respect to 
N. 



evaluate the curl and the divergence of n. The curl of n 
is calculated as 



Y, = a,Y, Yy = 9,9jY = F^.Yfc + L,,N, 
g,j=Yi-Yj g*-' EE (gy)"\ gEEdetgy, 



V X n = 



L, 



Y,, • N, 






(7) 



V5 



{gikdjHk + guUkT'-j^.) N 
- {2nkLjk + dj cos 6) Yi] , 



(13) 



{LvY 



L = det L, 



where eijk is the antisymmetric Levi-Civita tensor. Us- 
ing Eqs. (ITU)) and ([T^ we obtain the twist of the director. 



n • (V X n) = -^ [cos 6* {gikdjUk + garikT^^.) 

- {2nkniguLjk + gunidj cos6)] 



The indices i,j,k — 1,2 (repeated indices being summed 
over), di = 8^, and dj\[ denotes the partial derivative 
in the normal direction N. The completely symmetric 
tensors gij and Lij are the first and second fundamental 
forms of the surface, respectively. The Christoffel sym- 
bols F^^- are defined by the relation F^^- = g'""^Yij ■ Y,„, 

and obey the symmetry property F^- = F^;. The unit The divergence of the director is given by 
normal vector of the surface is given by 

Y1XY2 



V • n =: djfij + niT-'-i — cos OLjig'^ 



N 



V9 



The Gaussian curvature Kq, is given by: 



Kr,^-, 



9 



(8) 



(9) 



= 9jnj 



niT'^i , 



(14) 



(15) 



where we have used the fact that H ~ ^g^^ Lij = for a 
twisted ribbon Using Eqs. ©, dH]), (O, (HI]), the free 
energy Eq. ([T]) is given by. 



and the surface area element is given by dA = y/gduidu2. 
The director field n can be expressed in the local basis 
{Yi,Y2,N}as: 



F 



K 



{djn,+niT]iY 



n ~ riiYi + cos 6*1^ , 



(10) 



where the director components n^, i = 1,2 and the tilt 
angle are functions of the surface coordinates ui, U2 (see 
Fig. 1). 

The unit length of the director field yields the con- 
straint: 



-2?-^ [cos 9 {gikdjUk + gunkV^jk) 
v9 
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sin^ 6* = . 



(11) 



-(- — / sin 0y/^duidu2 + k / —y^duidu2 
+ ^ (p dl. 



The three-dimensional gradient operator is given by 



V - g'^Y.dj + N9 



'N ■ 



(12) 



where g^^Yidj describes the components of the gradient 
operator on the surface. To evaluate Eq. (jH]), we need to 



As noted in Ref. [T^, the term proportional to q in the 
equation above is identical in form (up to a line integral 
that can be absorbed into the edge energy term propor- 
tional to 7) to the Helfrich-Prost term [IJl first intro- 
duced to model chirality in lipid membranes. 



B. Sm-j4* twisted ribbon with Gaussian curvature 

term 

The position vector of a twisted ribbon (Fig. 2) of ra- 
dius R is given by 



Y — {p cos (p, p sin 0, 60} , 



(17) 



where \p\ < R and Q < (j) < 2mr, and n is the winding 
number. The pitch of the twisted ribbon is given by 
27r|6|, and the sign of 6 determines the handedness of the 
ribbon. Here we focus on a right-handed ribbon with 
6 > without any loss of generality. 

Defining ui = p and U2 = (j>, the tangent vectors and 
the normal vector of the surface are given by 



Yi = dpY ~ {cos 0, sin 0, 0} , 
Y2 — d^Y — { — p sin 0, p cos 0, 6} 

N = {6 sin (j), —6 cos 0, pi . 

Using Eq. dH]) in Eq. ^ we find: 
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(19) 



As we will ultimately carry out a mean-field analysis of 
the free energy, we assume that the director field exhibits 
azimuthal symmetry. Furthermore we assume that the 
director is tilted parallel to the edge of the ribbon. The 
director field n is then given by 



sin 6'(p) 
n= <jO, — ,cos6'(p) 



(20) 



where we have used the unit length constraint, Eq. ([TT|). 
Note that the director field of Eq. (PUJ) has zero splay; 
this can be verified explicitly using Eq. (ITSl) . 

Within the single Frank constant approximation we 
have one penetration length given by At = (■^)^'2. The 
thickness of the membrane is included in the elastic mod- 
ufi K and C. Inserting Eqs. ([H]), dH]), and (gHl) into 
Eq. (J16p we find the following expression of the free en- 




FIG. 2: A twisted ribbon with radius R showing the director 
field determined by the solution of the Euler-Lagrange equa- 
tion, Eq. (|22[) . The figure shows a length of the ribbon equal 
to 0.75 times the pitch. 
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The Euler-Lagrange (EL) equation for 9 is given by: 
d26l(p) ^ ^^d9{p) _ /'Zl^ , g 

g 



dp^ 
A\ pq 



dp 
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sin2 6l(p) = . 



The boundary conditions are given by: 
9{Q) = 

( d9{p) 



2{R'' + V) 



\ dp 



psm29{R) 



p=R 



+4bsm^9{R) =0. 



(22) 



(23) 



(24) 



Eq. (|23l) arises from the relation 9{p) — —9{—p) on a 
twisted ribbon, and is also observed experimentally, while 
Eq. ([M|) is the free boundary condition, ^^ = 0, where 

9' = ^. Physically, this boundary condition ensures 
that the director torque vanishes at the edge of the rib- 
bon. The director field given by the solution of Eq. (|22p . 
subject to the above boundary conditions, is shown in 
Fig. 2. 

The EL equation and boundary conditions derived in 
Ref. 4] for a semi-infinite smectic layer (Fig. 3) with a 



',;v\;v\:n;n:\;\:\;\:\;\:\\'\\\\\:\;\\\:'' 



\sssimw-^'' 



>>;'jW^ 



^?ivmi^ 



ym22^^ 



m^- 



FIG. 3: Director field near the edge of a semi-infinite smectic 
layer, the limit of a twisted ribbon with & — >■ oo. 



straight boundary can be recovered from the above re- 
sults for the twisted ribbon by taking the limit 5 -> oo. 
On the other hand, the limit 6 — >■ with fc = and R 
constant results in a finite disk of radius R. This limit 
leads to the EL equation and boundary conditions for a 
finite disk given in Ref. [8] . 



C. Sm-A* twisted ribbon with saddle— splay term 

We now consider a modification of the above free en- 
ergy of a ribbon replacing the Gaussian curvature term 
in the Helfrich energy Eq. ([2|) by the Frank saddle-splay 
energy density, i.e, we replace the surface normal N by 
the director field n in Eq. ([3]) and obtain the saddle-splay 
free energy density: 



h 



Ko 



-V- (n(V • n) - (n- V)n) , 



(25) 



where ^24 is the saddle-splay modulus. Two major sim- 
plifications help us evaluate Eq. ((25)) : (1) the director 
given in Eq. (|20|) has no radial component; (2) due to 
our choice of the coordinate system, metric tensor gij of 
the surface given in Eq. p^ is diagonal. Hence, Eq. P^ 
simplifies to 



V = g"Yi9i + g22Y2a2 + Na. 



N ■ 



(26) 



Calculating the first term on the right-hand side of 
Eq. ([^5)) apart from the factor ^^, we obtain 



V • (n(V -11)) = -AH^ cos^ 9 + 2Kg cos^ 
while the second term of Eq. (P5|) yields 



(27) 



V • ((n • V)n) =di {{n2fTl2) - 2<9i (cos6'n2L2i) 
+ {n2fTl2d2Yi -2cos9n2L2iTl^ 
- ("-2^21)^ , 



(28) 



with the director component 712 given by Eq. (l20l) . Fi- 
nally, the saddle-splay free energy density, Eq. ([25]) is 
given by: 
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(29) 



-b ( 2gcos26'- p sin 20 



In the limit of zero tilt angle (i.e., n = N), Eq. (|29|) 
reduces to the term proportional to k in Eq. (f2T|) if we 
replace /'ir24 by k. The EL equation including the saddle- 
splay term in place of the Gaussian curvature is given by: 
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(30) 



with the boundary condition at the edge p = R: 
- 2 {qR^ + b{bq - 1)) - 26(^ + 1) cos 261 
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+ {^ + l)Rsm20 + 2{R' + b')^ 
K op 



(31) 



= 0. 



p=R 



III. RESULTS 

We now explore the results of solving the EL equa- 
tions and their respective boundary conditions for the 
two models of the twisted ribbon: with the Gaussian cur- 
vature term (Eqs. ([22]) -([M l) ) and the saddle-splay energy 
(Eqs. (15(11) -([5T |) V We ignore any change in the membrane 
thickness due to the tilt of the director and assume a com- 
mon value of the Gaussian modulus k and the saddle- 
splay modulus ^24 in comparing the two models. For 
a given value of these parameters and 7, the edge en- 
ergy, we determine, using a downhill-simplex algorithm, 
the values of the radius R and pitch b which minimize 
the free energy of the ribbon. We then compare the 
free energy per unit area of the twisted ribbon with the 
corresponding energy density of a semi-infinite Sm-A* 
membrane [4j, where as in the case of the ribbon, the 
director is tilted parallel to the edge of the membrane 
(Fig. 3). We determine the phase boundary between the 
two structures using a bisection algorithm. The area A^ 
of a twisted ribbon is given by 



Ar= I ^fgd-pd^p = 2ti I yjp^ + b'^dp 
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FIG. 4: Phase diagram for semi-infinite membranes and 
twisted ribbons witli q\t — 0.71. The results are shown both 
in dimensionless units (left and bottom axes), K24 = -^ 
and 7 = 2_i^ and physical units (top and right axes) ob- 
tained using the value of At measured in Q| and the value of 
K measured in [6|. The solid curves denote first-order phase 
transition lines for the Gaussian curvature (black) and saddle- 
splay (gray) models. The saddle-splay modulus K24 and the 
Gaussian modulus k are taken to be equal. The dashed line 
is not a phase boundary; rather it denotes the second-order 
phase transition between finite-sized disks and semi-infinite 
membranes obtained in a theoretical model [g which did not 
consider twisted ribbons. The twisted ribbon is of lower en- 
ergy than the finite-sized disks. 



The energy per unit area of the semi-infinite membrane 
is given simply by q'^/2 as director tilting occurs only 
in a small region near the edge, and thus in Eq. ([5]) the 
contributions of twist and bend deformations to the mean 
free energy density are negligible. 

We have considered both positive and negative values 
of the moduli k and K24- We find that negative values of 
these quantities lead to stable ribbons with ratios of pitch 
to radius of order 35 or greater, significantly larger than 
what is measured experimentally [0]. On the other hand, 
positive values of the moduli lead to stable ribbons with 
ratios of pitch to radius of order 5, in good agreement 
with experiment. Positive values of these moduli are not 
typically measured in lipid monolayers or bilayers [18| : 
however, there is no reason to exclude this possibility on 
physical grounds and the fd system may be very different 
from systems composed of amphiphilic molecules. Posi- 
tive moduli lead to saddle-splay and Gaussian curvature 
energies which are negative for a twisted ribbon. Thus, 
to stabilize the free energy we have added a higher-order 
curvature term proportional to T^Q [19iji2C|]- The coupling 
of this term to the free energy can be chosen as small as 
0.01 in dimensionless units to guarantee stability for the 
range of K24 and k values we have studied. 

Our main results are displayed in Figs. 4-6, where we 
present the phase diagram for semi-infinite membranes 
and twisted ribbons, and values of the tilt angle, pitch 
and width of the ribbon as functions of dimensionless 



variables ii'24 



-^ and 



7 



K 



-, as well as physical 



units obtained by using the value of K measured in [6| 
and At measured in [8|. Our results were obtained for 
qXf = 0.71, corresponding to the values of these parame- 
ters measured experimentally in Ref. [8[. We find a first 
order phase transition between ribbons and semi-infinite 
membranes (Fig. 4) as the edge energy 7 (proportional to 
the depletant concentration) is varied. Furthermore, we 
see that even though there is a finite tilt of the molecules 
throughout the layer (see Fig. 5c) which should in prin- 
ciple distinguish between the saddle-splay and Gaussian 
curvature models, there is no significant difference be- 
tween the calculated phase boundaries for the two mod- 
els. The dashed line in Fig. 4 is not a phase bound- 
ary; rather it is the second-order phase transition line 
between semi-infinite membranes and finite-sized disks 
calculated in a model |8| which excluded the possibility 
of twisted ribbons. The present calculation demonstrates 
that ribbons are of lower free energy than the finite-sized 
disks and the first-order phase transition from semi- 
infinite membranes to ribbons preempts the transition to 
finite-sized disks as 7 is reduced, in agreement with ex- 
perimental observations [9[ . An examination of the var- 
ious contributions to the free energy indicates that it is 
the reduction in the chiral energy density —qKn- (V x n) 
that drives the transition to the ribbon state. If q, the 
measure of chirality, is set equal to zero, ribbons are never 
the lowest free energy structure. 

The upward slope of the phase boundary in Fig. 4 in- 
dicates that as 7^24 increases the ribbon is energetically 
more favorable than the semi-infinite membrane. This 
tendency arises because increasing i^24 makes the saddle- 
splay energy (or equivalently the Gaussian curvature en- 
ergy) more negative, thus making the ribbon more favor- 
able. This effect outweighs two factors making ribbons 
less favorable as K24 increases, namely, the bend energy 
increases and there is less twist penetration at the edge. 
The increase in bend energy can be inferred from Fig. 6c 
where the ribbon pitch b is plotted as a function of ^24 
for fixed 7. Recalling (see the end of Sec. IIIBp that the 
limit 6 — >■ 00 corresponds to a semi-infinite membrane 
(where the bend deformation vanishes) while the oppo- 
site limit 6 — > corresponds to a disk of radius it! where 
the director field exhibits bend f8| , Fig. 6c indicates that 
the bend energy increases with increasing i^24- The de- 
crease in twist penetrations as -ftr24 increases can be seen 
from Fig. 5a where the tilt angle Oq = 0{R) at the edge 
of the ribbon is plotted as a function of A'24. 

While Fig. 4 is in qualitative agreement with exper- 
iment, namely, that a first order transition is observed 
between semi-infinite membranes and twisted ribbons, 
the predicted line tension at the phase boundary, 7 ^ 
10 kBT/fj,m is an order of magnitude lower than the value 
measured in experiments [9|. This discrepancy persists 
even when one accounts for the fact that the experiments 
measure an effective line tension, i.e., including the net 
reduction of the free energy near the edge due to twist 
penetration. 

The tilt angle 6*0 = d{R) at the edge of the ribbon for 
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FIG. 5: (a) The value of the tilt angle 60 = 9{R) at the edge of the ribbon as a function of K24, for 7 — 0.269, corresponding to 
a point slightly below the first-order phase boundary shown in Fig. 4. The black and gray curves correspond to the Gaussian 
curvature and saddle-splay models respectively. The values of 60 shown are all less than the value found for the semi-infinite 
membrane ,4]; (b) 60 as a function of 7 for K24, ~ 0.15, as the phase boundary is approached from below. The vertical dashed 
lines denote from left to right the first-order phase transition from twisted ribbons to semi-infinite membranes for the Gaussian 
curvature and saddle splay models. To the right of these lines the ribbon is metastable; this is indicated by the dashed nature 
of the do curves, (c) The tilt angle 9{p) for the saddle-splay model with K24, = 0.15 as a function of p in units of At (lower 
ajcis) and physical units (top axis) corresponding from top to bottom to the parameters: 7 = 0.25, {R — 2.566At, b = 2.877At); 
7 = 0.26, (R = 2.739At, 6 = 3.381At,); 7 = 0.27, {R = 2.965At, b = 3.899At); and 7 = 0.28187, {R = 3.385At, b = 4.707At,). 
The lowest curve (7 = 0.28187) corresponds to the phase boundary. 



both the Gaussian curvature and saddle-splay models is 
shown as a function of both K24 and the edge energy 7 
in Figs. 5a and 5b respectively. From Fig. 5a we see that 
as with the phase boundary shown in Fig. 4 there is neg- 
ligible difference between the results for the two models. 
The vertical dashed lines in Fig. 5b correspond to the 
first-order phase transition boundaries shown in Fig. 4 
for the saddle-splay and Gaussian curvature models. To 
the right of these lines the ribbon is metastable which we 
indicate by the dashed plotting of Oq. For semi-infinite 
membranes 3 610 = — arcsin((7At) = 0.79 for qXt = 0.71, 
independent of 7. This value exceeds Oq for the twisted 
ribbon in both the saddle-splay and Gaussian curvature 
models for all values of 7 and K24, in agreement with ex- 
periment Q . In Fig. 5c we show examples of the tilt angle 
profile 6{p) for the saddle-splay model with K24 = 0.15 
for a range of 7 slightly below and on the first-order 



phase boundary of Fig. 4. 

The pitch 27r6 and width 2R of the ribbon are shown as 
functions of 7 and K24 in Figs. 6a-d. In Figs. 6a and 6b 
we see that there is some difference between the Gaussian 
curvature and saddle-splay models for the dependence of 
the pitch on 7; we have no simple physical reason for this 
difference. Fig. 6e shows the geodesic curvature Kg [ij| 
of the edge of the ribbon: 



R 



R^ + b^ 



(33) 



as a function of 7 in the neighborhood of the first-order 
transition. For both the Gaussian curvature and saddle- 
splay models the geodesic curvature decreases as the 
phase transition is approached from below, while Fig. 5b 
shows that the ^o increases. This result is consistent with 
experimental retardance data which shows that lower 
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functions of 7 at K2i = 0.15. Black and gray curves correspond to the Gaussian curvature and saddle-splay models respectively. 
The gray vertical dashed lines denote from left to right the first-order phase transition from ribbons to semi-infinite membranes 
for the Gaussian curvature and saddle-splay models, respectively. To the right of these lines the ribbon is metastable; this is 
indicated by the dashed nature of the curves. In (c) and (d) respectively the pitch and width are shown as functions of K24, 
close to the phase boundary at 7 = 0.269. Again, black and gray curves correspond to the Gaussian curvature and saddle-splay 
models respectively. The geodesic curvature, Eq. ^i'i\ . as a function of 7 is shown in (e). 



curvature of the edge in general leads to a higher tilt 
angle at the edge. 



IV. CONCLUSION 

In this paper, we have developed a theoretical model 
for the Sm-^* twisted ribbons observed in assemblies of 
fd viruses condensed by depletion forces. The depletion 



interaction is modeled by an edge energy (line tension), 
assumed to be proportional to the depletant polymer in 
solution. We have considered two variants of this model, 
one with the conventional Helfrich Gaussian curvature 
term, and a second with saddle-splay energy. Both mod- 
els yield a first-order phase transition between ribbons 
and semi-infinite fiat membranes as the edge energy is 
varied. The phase transition line and tilt angle profile 
are found to be nearly identical for the two models; the 
pitch of the ribbon, however, does show some differences. 
Our model yields good qualitative agreement with exper- 
imental observations, namely, (1) the existence of a first- 
order phase transition between ribbons and semi-infinite 
membranes, preempting a transition to finite-sized disks; 
(2) the dependence of the tilt angle at the edge of the rib- 
bon on the curvature of the edge; (3) the decrease in the 
tilt angle at the edge of the membrane as the system un- 
dergoes the first-order transition from the semi-infinite 
membrane to the ribbon; (4) the order of magnitude of 
the pitch to width ratio of the ribbon. While the phase 



diagram is in qualitative agreement with experiment the 
predicted line tension on the phase boundary is an order 
of magnitude lower than that measured experimentally. 
Also, as we noted in Sec. IIIIl obtaining the correct or- 
der of magnitude of the pitch to width ratio requires 
a positive value of the Gaussian curvature or saddle- 
splay modulus which favors negative Gaussian curvature. 
These moduli have not been directly measured in the fd 
system nor have microscopic models been developed for 
nonamphiphilic molecules. 



Acknowledgments 

We thank E. Barry, Z. Dogic, T. Gibaud, and M. Za- 
khary for useful discussions. This work was supported 
by the NSF through MRSEC 0820492 and by Brandeis 
University. 



[1] P. G. De Gennes and J. Prost, The Physics of Liquid 

Crystals (Oxford University Press, 1993). [12] 

[2] P. G. De Gennes, Solid State Commun. 10, 753 (1972). [13] 

[3] See, e.g., R. D. Kamien and J. V. Selinger, J. Phys: Cond. [14] 

Mat. 13, Rl (2001) 
[4] E. Barry, Z. Dogic, R. B. Meyer, R. A. Pelcovits, and R. [15] 

Oldenburg, J. Phys. Chem B 113, 3910 (2009). 
[5] E. Barry, D. Beller, and Z. Dogic, Soft Matter 5, 2563 [16] 

(2009). 
[6] Z. Dogic and S. Fraden, Langmuir 16, 7820 (2000). [17] 

[7] Z. Dogic and S. Fraden, Curr. Opin. Colloid Interface Sci. [18] 

11, 47 (2006). [19] 

[8] R. A. Pelcovits and R. B. Meyer, Liq. Cryst. 36, 1157 

(2009). [20] 

[9] T. Gibaud and Z. Dogic, private communication. 
[10] W. Helfrich, Naturforsch. 28, 693 (1973). 
[11] W. Helfrich and H. J. Deuling, J. Phys. (Paris) CoUoq. 



36, C1327 (1975). 

Z. C. Tu and U. Seifert, Phys. Rev. E 76, 031603 (2007). 

See, e.g., R. D. Kamien, Rev. Mod. Phys. 74, 953 (2002). 

M. Kleman and O. D. Lavrentovich, Soft Matter Physics: 

An Introduction (Springer Verlag, 2003). 

Ou-Yang Zhong-can and L. Ji-xing, Phys. Rev. Lett. 65, 

1679 (1990), Phys. Rev. A 43, 6826 (1991). 

Ou-Yang Zhong-can and W. Helfrich, Phys. Rev. A 39, 

5280 (1989). 

W. Helfrich and J. Prost, Phys. Rev. A 38, 3065 (1988). 

D. Marsh, Chem. and Phys. of Lipids, 144, 146 (2006). 

T. Shemesh, A. Luin, V. Malhotra, K. N. J. Boerger and 

M. M. Kozlov, Biophys. J. 85, 3813 (2003). 

M. D. Mitov, Comptes Rendus Acad. Bulgares Sci 31, 

513 (1978). 



